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We propose a mechanism to make gravitational waves (GWs) visible in the electromagnetic domain. Grav- 
itational waves that propagate through a strongly magnetized plasma interact with the plasma through its 
anisotropic stress-energy tensor and excite magnetohydrodynamic (MHD) wave modes. In catastrophic events 
such as the merger of a double neutron star binary, a large fraction of the total binding energy of the system is 
released in the form of GWs observable by LIGO, and the amount of energy transferred to the MHD waves is 
substantial. These modes, however, are excited at the same frequency as the GW and are not directly observ¬ 
able. In this paper we investigate radiation processes that operate in the presence of the gravitationally excited 
MHD waves and radiate in the radio regime accessible to LOFAR. We present order of magnitude estimates for 
the spectral flux of a merger detectable by a LOFAR. 


1. INTRODUCTION 

Gravity is the weakest of all the fundamental forces 
and, not surprisingly, gravitational waves have proven 
to be the most elusive signals in astronomy. In fact, 
now that almost every electromagnetic frequency is 
observable and even cosmic ray and neutrino astron¬ 
omy are becoming a reality, the first direct detection 
of GWs still has to be made. With the development 
of advanced resonant bar and in particular interfer¬ 
ometric GW detectors, this will most likely happen 
within the next one or two decades. The same prop¬ 
erty of GWs that makes them so hard to detect is their 
main advantage as a diagnostic of the Universe: be¬ 
cause GWs hardly interact with anything the entire 
Universe is transparent to the m. The planned GW 
detector in space, LISA IltsaI expects to detect the 
signal from merging su permas sive black holes out to 
any distance and LIGO TltgoI which is presently un¬ 
der construction could probe for instance the central 
engine of gamma-ray bursts (GRBs) which is opaque 
to electromagnetic radiation. 

In general, electromagnetic and GW observations 
are very complimentary and an electromagnetic coun¬ 
terpart of any GW signal and vice versa is invaluable. 
As an example: with electromagnetic observations it 
is very difficult to determine the inclination of a bi¬ 
nary system, whereas this could be easily measured 
from the polarization of a GW signal ISchut d jl990j . 
On the other hand, the amplitude of a GW measured 
on earth does not allow a unique determination of the 
distance and the intrinsic strain at the source (which 
can be measured only from the chirp of the system). 
When the distance can be determined electromagnet- 
ically, this uncertainty can readily be removed. 

Rather than studying different types of GW sources 
and trying to find accompanying electromagnetic 
processes, we investigate whether the GW them¬ 


selves produce an electromagnetic signal when they 
propagate through the magnetohydrodynamic (MHD) 
plasma in which many GW sources are embedded. In 
particular one can think of rapidly spinning NS with 
small (< 10” 5 ) asymmetry, low-mass X-ray binaries, 
neutron stars with an r-mode instability, asymmetric 
supernova core collapse and bounce (GRB with after¬ 
glow), newly born neutron stars that boil and oscillate 
(magnetars) and coalescing compact binaries (short 
GRB candidates). Gravitational waves only interact 
to linear order with mediums that have an anisotropic 
stress-energy tensor. Interaction is therefore possible 
with a magnetized MHD plasma and not, for instance, 
with dust or an isotropic hydrodynamic ideal fluid. 

We ha ve fo und in pr evious work 

(see _[Moort gat andKuiipersI l2003| . 

iMoortgat and Kubpe^T2004j and references therein) 
that x and + polarized GW couple to Alfven and fast 
magneto-acoustic MHD modes, respectively. Here 
the x polarization refers to GWs that are polarized 
at an angle to a uniform magnetic field. In the limit 
of a strongly magnetized tenuous plasma these MHD 
waves can interact coherently with a GW and grow 
linearly with distance up to a maximum value. The 
frequency of the excited MHD modes is, however, 
the same as that of the GW and therefore lower than 
the interstellar plasma frequency and not observable 
directly. 

In this paper we investigate subsequent inverse 
Compton (IC) scattering of the low frequency 
(kHz) MHD wave on the relativistic electrons (and 
positrons) in the wind, either the bulk secondary par¬ 
ticles or the even more relativistic primary particles. 
Inverse Compton scattering on electrons and positrons 
with a Lorentz factor of 7 > 100 boosts the f requency 
of the radiation to the regime of the LOFAR iLOFARl 
low frequency radio array currently under construc¬ 
tion in the Netherlands. Finally, we mention the possi- 
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bility of coherent radiation when the MHD mode acts 
as an undulator for the relativist ic electrons similar 
to the free electron laser process iFung and Kuiipersl 
200 1 1 . 

The set-up of this paper is as follows: Sec¬ 
tion O sum marizes some of the main con¬ 
clusion^ of_ Moor tgat and Kuiipersl j2003j and 
iMoortgat and Kuiipersl [2004j . Section 0 elaborates 
on the specific situation of GWs propagating through 
the magnetosphere and wind of a double neutron star 
(neutron star) binary close to merging. In Sections 0] 
- Section^we derive a rough estimate for the spectral 
flux on earth as a result of the inverse Compton 
scattering, taking into account possible damping 
mechanisms and propagation effects and we discuss 
some issues that are subject of future research such 
as the coherent radiation mechanism. We end with 
conclusions in Sectional 


B x °+ Bj 1 (z) 


-2 0 2 4 



Figure 1: Magnetic field in the magneto-acousic mode. 
The perturbation is exaggerated to emphasize the overall 
behavior (arbitrary units). 


2. MHD MODES EXCITED BY GW 


The only oscillating components 5T ab of the stress- 
energy tensor in the rest frame of a perfect mag¬ 
netofluid that couple to the GW and can not be re¬ 
moved by a gauge transformation depend on the mag¬ 
netic field. Explicitly, the Einstein field equations re¬ 
duce to evolution equations for the GW: 

nh+{z,t) = 4 B°5B x (z,t), (1) 

nh x {z,t) = 4 B°JB y (z,t), (2) 


where the ambient magnetic field is chosen to lie in the 
x-z plane: B = B° + SB and B° = B °(sin 0, 0, cos 0). 

Approximating the GW as a driver, the evolution 
equations for the magnetic field can be solved and look 
like: 


SB X oc - h + B °, 


6B y oc -h x B° x 


( 3 ) 


These results are reminiscent of the spatial deviations 
of test masses in interferometers such as LIGO (8x = 
\{h+x o + h x y 0 ) and Sy = \(h x x 0 - h+y 0 )). 

The solution for 5B X is illustrated in Figure Q 
and corresponds to a compressional fast magnetosonic 
wave (MSW) with both electromagnetic and gas prop¬ 
erties. Coherent interaction with the GW is possi¬ 
ble when the phase velocity of the MSW approaches 
that of the GW. In a Poynting flux dominated plasma 
where the Alfven velocity ua is relativistic and much 
larger than the sound velocity, this limit is satisfied 
and the perturbations are allowed to grow linearly 
with distance: 


SB x (z,t) ~ -^-B 0 sin0 coz 7s 

4-y2 


j(z-t) 


( 4 ) 



Figure 2: Magnetic field in the Alfven mode. 


The second expression in Eq.|3] corresponds to non- 
compressional shear Alfven waves, illustrated in Fig¬ 
ure |2J The condition for coherent interaction with 
the GW is more stringent because its phase velocity 
U A|| = ua cos 6 has to approach the velocity of light, 
but at the same time its amplitude is oc B x oc u\ sin0. 
Therefore, in the case of coherent interaction the am¬ 
plitude of the Alfven waves is suppressed by a small 
factor 0 <C 1. 


5B y {z,t) ~ ^B°0 toz S[e i “ (z_t )] + O[0 2 ]. 


( 5 ) 


The phase velocity of the slow MSW is always much 
smaller than the fast mode, so it can never interact 
coherently with the GW. 


In fact the polarization of the GWs emitted by a 
binary merger is linear only for an edge-on and circular 
for a face-on binary with the same sense as the binary 
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Figure 3: Polarization of Alfven (top) and fast 
magnetosonic modes (bottom). 


motion on the sky. If the GW is for instance right 
circularly polarized, simply replace h+ —> /ir/v 7 2 . 

As was mentioned in the previous section, the GW 
only interacts directly with the magnetic field. In 
particular the plasma motion in a GW is generally 
non-compressional. However, in a perfectly conduct¬ 
ing plasma the particles are strongly coupled to the 
magnetic field lines and Maxwell’s equations couple 
to the matter conservation laws through the current 
density. Consequently, perturbations in pressure, den¬ 
sity, magnetic field gradients, currents and a drift ve¬ 
locity are also excited in the MSW, whereas the non- 
compressional Alfven waves cause a divergence of the 
electric field and a corresponding charge density fluc¬ 
tuation. The orientation of these different components 
of the MHD modes, i.e. the polarizations, are shown 
in Figure [3J 


3. GENERAL PICTURE 

We consider a merger of two neutron stars where 
each neutron star has its own magnetic field that falls 
off as a dipole close to the surface: B(r) = B* (i?*/r) 3 . 
As the merger coalesces the orbital frequency increases 
and dominates over any other (rotational) motion. 
Therefore, we assume that at the end of the spiral- 
in phase the orbital rotation of the binary (with 
fib ~ 10 3 rad/s) determines the light-cylinder ra¬ 
dius: R\ c = c/flb — 320 km. At the light cylinder 
plasma would require superluminal velocities to coro¬ 
tate. A steady charge density cannot be maintained 



Figure 4: Magnetic field configuration of a double 
neutron star binary merger. 


at the Goldreich-Julian density, uq .t = 2eo VL h -B/e, 
everywhere. As a result, a strong electric field devel¬ 
ops along the open magnetic field lines above the po¬ 
lar cap and charged ‘primary’ particles are extracted 
from the surface with a density n p — 71gj* and ac¬ 
celerated to high Lorentz factors (a typical number is 
7 P ~ 10'). The available potential jump is propor¬ 
tional to oc flb-B* in our case and can be much 
larger than for single pulsars since the orbital fre¬ 
quency becomes large for neutron stars of arbitrary 
field strengths. 

Processes such as curvature radiation and inverse 
Compton emission then result in a cascade of ‘sec¬ 
ondary’ pairs with a particle number density n s = 
2 Mn p , where M is called the multiplicity. Due to 
energy conservation n p y p = ro s 7s, so the Lorentz fac¬ 
tor of the secondary particles is 7 S = 7 P /( 2 M) ~ 100 
for M ~ 10 5 . The secondary plasma, probably in¬ 
terspersed with beams of primary particles, flows out 
as a relativistic wind along the open magnetic field 
lines which develop into spirals further out since the 
toroidal component of the field soon dominates the 
poloidal component. The wind remains force-free up 
to a large distance. 

As the binary looses angular momentum through 
gravitational radiation it spirals in and eventually the 
two neutron stars merge and form a black hole. In 
the process, a large fraction of the binding energy of 
the system is released in the form of GWs that propa¬ 
gate through the surrounding plasma and excite MHD 
waves. In the relativistic wind the dominant mode will 
be the fast magnetosonic wave that propagates with a 
highly relativistic Alfven phase velocity u\ | c, over¬ 
taking the particles further out (that are themselves 
relativistic). 

In the next section we will discuss how the energy 
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transferred to the plamsa by the GW can be released 
in the form of inverse Compton radiation. 


4. INVERSE COMPTON RADIATION 


The relativistic flow in the jet cannot remain force 
free after the GW pulse has travelled through the jet 
and excited the MHD waves which - by definition - 
impose forces and accelerations everywhere on the lo¬ 
cal plasma. As a result the jet is expected to radiate. 
For simplicity we assume that the particles are charac¬ 
terized by an outflow of secondary particles at multi¬ 
plicity M and a random motion with Lorentz factor 7 S 
superimposed on a systematic outflow with a similar 
Lorentz factor. Of course, there will also be the rem¬ 
nant of the primary beam with Lorentz factor but 
these particles will not radiate at radio frequencies. 
We note however that the secondary particles will (see 
Eq. [HI below). Actually, the radiation from the sec¬ 
ondary particles will be partly synchrotron in the mag¬ 
netic field of the MHD waves - largely at frequencies 
above the LOFAR band - and partly inverse Comp¬ 
ton radiation. Although the inverse Compton process 
is fixed by the frequency of the GW and the Lorentz 
factor of the particles - and has therefore the same fre¬ 
quency and in the radio domain all over the jet - this 
radiation can only escape sufficiently far out where 
the relativistic plasma frequency has fallen below the 
frequency of emission oj s > u p /\J < y 3 > ~ u> p /y/"U 
where u> p is the non-relativistic plasma frequency of 
secondary particles in the observer frame and the av¬ 
erage is over the energy distri b ution of the secondary 
particles iMelrose and Gedalinl Il999l . 

Inverse Compton scattering is a natural radiation 
process because the GW - and therefore also the MHD 
wave - propagate outwards radially, whereas the wind 
is collimated by the background magnetic field and 
has an essentially cylindrical geometry. Consequently, 
the MHD wave overtakes the wind at an angle, and 
in fact the relativistic electrons will see it approaching 
almost head-on in their rest frame (Figure El). In fact, 
since we are studying the regime of a strongly magne¬ 
tized tenuous plasma, the MHD wave behaves similar 
to a vacuum electromagnetic wave propagating at the 
speed of light. We will therefore treat the MHD wave 
simply as incident photons with to = w m hd = w gw in a 
first estimate. 

The frequency of the scattered radiation can most 
readily be derived from the resonance condition: 


^rn hr] COS O gw — 0J S k s /3 COS (X s 



where the last equality is for the photons scattered in 
the forward direction (a s = 0). The frequency of the 
scattered radiation as a function of the angle between 
the incident MHD wave and the bulk velocity in the 
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Figure 5: Inverse Compton process in plasma (left) and 
lab (right) frame. 


wind is plotted in Figure El Already for relatively 
small angles a gw does the boosted frequency fall into 
the LOFAR regime, which is 10 - 220 MHz. 



Figure 6: Scattering of the radial MHD waves on the 
cylindrical flow in the wind (in this illustration the angle 
ogw is compressed). 


The p ower dens i ty in inverse Compton radiation is 
given by iRvbicki and LightmarJ Il979l 

Pic ~ K7 s 2 /3 2 p T eEMn e , (7) 

where n e is the density of scattering electrons (in fact, 
the scattered radiation from the electrons and the 
positrons adds constructively). f3 is the velocity of the 
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electrons and 7 S the corresponding relativistic Lorentz 
factor, £em is the energy density in the incident MHD 
waves and k is a geometrical factor of order unity, 
which depends on the angular distribution of the in¬ 
cident radiation (k = 4/3 for an isotropic radiation 
field). Note that the Thomson cross-section, or, is 
only a rough estimate for the scattering of a MHD 
wave on relativistic charged particles. 

The energy density in the MHD mode (with SE 2 + 
SB 2 ~ 25B 2 \pdv 2 ) is iMoortgat and Kuiipersl 
12004 1: 


£em 


\_SBl 

Oo 


' (B °) 2 1 (k mhd z) 2 h 2 + 

. 2^ 0 J 2 7 4 


( 8 ) 


For a GW burst of duration At < 0.1 s, 

(fcmhdZmax) 2 fr+ = Lo 2 h\At 2 = 0[ 1] and the excited 
MHD wave amplitudes can be quite significant for a 
mildly relativistic wind. Note also that in the analyt¬ 
ical calculations we have assumed a uniform constant 
background magnetic field. In the numerical estimates 
in the next section, we will take account of the 1/r de¬ 
pendencies in a WKB approach so that the MHD wave 
amplitude falls of as 1 /r within the interaction region 
just like the background field. 


4.1. Numerical estimates 

Radiation from the entire jet will be emitted as soon 
as the GW pulse has passed through and the MHD 
waves are excited. The burst of electromagnetic radi¬ 
ation will be compressed to a duration At ss L/(c2j 2 ) 
as is familiar from GRB models. However, even 
though the entire jet is radiating at the same inverse 
Compton frequency the jet expands and is therefore 
very inhomogeneous. In our estimate of the total 
power in inverse Compton radiation we will just ap¬ 
proximate the jet as being built up out of different 
parts each of uniform properties and extent equal to 
the local scale height L ss r. 

In the absence of collimation both the electron den¬ 
sity, the backgound magnetic energy density and the 
MHD wave field energy density fall off as 1 /r 2 and the 
simple estimate of Eq. 0 integrated over r 2 dr shows 
that the power falls of with 1/r. Radiation in the LO- 
FAR radio domain can only escape far away from the 
light cylinder where the plasma density has decreased 
sufficiently. Note however that a mild collimation of 
the jet with B oc 1 /y/r so that both CEM,n e oc 1/r 
already leads to the opposite result that the power 
increases (oc r). 

When the jet is collimated and using the values in 
Table Q] we find: 


p tot 
Me 


2,32 /' 2r eEM(-Rlc)^e(Rlc)^?i 2 c 
«7s P CT T / -3- 


r 2 dr 


4 • 10 30 W 



(9) 


Parameter 

Value 

B+ - 

o 10 12 T 

Hb ■" 

- 10 3 rad s _1 

7p " 

- to 7 

7s - 

- 10 2 

M - 

- 10 5 

Rlc = ^ r 

- 3 • 10 5 m 

ni — 2AfeoOh s l<, . 

- 3 ■ 10 23 m -3 

hi n r 

- 1(T 3 

<TT r 

- 66.5 • 10~ 30 m -2 


Table I Numerical values in Eq. 0 


Coherent interaction of the GW with the MHD is 
only possible as long as their wave number mismatch 
is small, (A k)z <C 1, from which we calculated a 
ma ximum i nteraction length scal e of i? max ~ 10 15 m 
in iMoortgat and Kuiipersl j2003l |. Most of the in¬ 
verse Compton radiation will be emitted in this last 
L = i?max with .Ff c ot ~ 10 4 ° W, which is approxi¬ 
mately 10 -5 of the power in GW! 

For a binary merger close enough for a LIGO detec¬ 
tion, i.e. in the Virgo cluster at D ~ 1 Gpc, this would 
yield a flux on earth of 2 ■ 10 6 Jy at 30 MHz assuming 
that the signal in reality is not monochromatic, but is 
incoherent with a bandwidth of Av ~ v = 30 MHz. 
The duration of this transient signal is relatively short: 
At = L/{2'jfc) « 3 min, but much longer than the 
cosmic ray events that LOFAR will also detect. 


5. DISCUSSION 

5.1. Propagation effects 

Whether this signal is sufficiently strong, depends 
on a number of complicating factors. Firstly, the 
galax y is quite bright at these low radio frequencies. 
From ICand Il979> we find that at 30 MHz the galac¬ 
tic background brightness temperature is 1.7 • 10 4 K. 
This is, however, still below LOFAR’s antenna tem¬ 
perature of ~ 2.3 • 10 4 K at the same frequency and 
well below our signal (both are approximately all sky 
tempe ratures an d can be therefore be compared di¬ 
rectly iPacholczvB I'l970j ). The galactic background 
temperature drops with frequency and it is even less 
of a problem at higher frequencies. 

Secondly, dispersion is important at low frequen¬ 
cies. Since the radio signal is propagating through 
both interstellar and intergalactic plasma rather that 
a vacuum, its dispersion relation becomes frequency 
dependent. Radiation at lower frequencies is delayed 
with respect to higher frequencies and a single pulse is 
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‘smeared out’ over a longer duration and with a cor¬ 
respondingly lower signal strength. This differential 
delay depends on the electron density, n e , integrated 
along the line of sight to a distance D, e.g. the elec¬ 
tron column density or dispersion measure DM in pc 
m -3 , so that the frequencies in Eq. 1111 are in MHz 
(and consequently also the differential dispersion is in 
s MHz- 1 ). 


dt 

dv 

DM 


-8 ■ 10“ 3 


DM 

IP' 



( 10 ) 

( 11 ) 


When looking at a neutron star binary merger the 
signal is dispersed both by the host galaxy and by 
the intergalactic medium. Since these binaries close 
to a merger are very old systems that have received 
a randomly oriented kick during the birth of each 
neutron star, we expect them to lie relatively high 
above the galactic plane where the galactic dispersion 
is low. The dominant dispersion then occurs during 
intergalactic propagation and results in a dispersion 
measure of DM~ 10' pc nT 3 for a source at a distance 
of 100 Mpc and an upper limit for the unknown inter¬ 
galactic electron density of n e ~ 0.1 m -3 . In terms of 
delays this means 


jl-3 s MHz -1 @ v = 30 MHz, 

dv , ’ ( 12 ) 

-0.08 s MHz" 1 @ v = 100 MHz. v ' 


For comparison: the dispersion measure in a typical 
galactic medium with n e ~ 3 • 10 4 m -3 over a distance 
of D ~ 3.5 • 10 3 pc is DM = 10 s pc m -3 , so if the 
merger occurs deeper inside the host galaxy, or the 
line of sight crosses our own galaxy close to the disc 
the dispersion is an order of magnitude larger. 

The effect of the dispersion on the signal will be 
to smear out the pulse within the observing band of 4 
MHz over say a minute. As a result the signal to noise 
will go down with the inverse square root of the in¬ 
crease in duration, but in view of the above estimates 
well within the sensitivity of LOFAR. 


5.2. Coherent radiation 

Here we mention the posibility of coherent radio 
emission from the jet. Once the fast mode MHD waves 
travel away from the jet axis into weaker magnetic 
fields their amplitude relative to the ambient magnetic 
field increases and may lead to the effective formation 
of a large amplitude wiggler or undulator. Just as in 
the Free Electron Laser the fast particle beam of sec¬ 
ondary particles emits coherent radiation, again at the 
frequency det ermined b y Eq. El b ut wit h an increased 
intensity (see iFung and Kuiiperi 12004 11. 


6. CONCLUSION 

We have shown that the pulse of GWs which is ex¬ 
cited in a merger sets up a pulse of fast mode outgoing 
MHD waves of the same frequency in the pre-existing 
jet structure of the binary. The interaction between 
fast (secondary) particles and these MHD waves cre¬ 
ates radio emission by the inverse Compton process 
which is expected to be observable with LOFAR. The 
radio emission pulse comes from the outer regions of 
the jet and is expected to be smeared out to longer du¬ 
ration by intergalactic dispersion but still detectable 
with LOFAR. The inverse Compton signal can serve 
as an indirect detection of GWs and be distinguished 
from other, more conventional radio emission from the 
merger, by its low frequency imprint. The best strat¬ 
egy to look out for these signals will be in GRBs and 
future observations of GW detectors such as LIGO. 
The coupling between GWs and MHD waves there¬ 
fore offers in principle an important new and indirect 
complementary detection technique of GWs. 
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